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ABSTRACT 

The  influence  of  Lewis  number  on  stability  and  oscillations 
for  stirred  tanks,  catalytic  wires  and  gauzes,  and  catalyst  par- 
ticles is  determined  quantitatively.  The  mathematical  structure 
is  discussed  and  illustrated  by  example.  Relaxation  oscillations 
are  found  to  be  the  natural  asymptote  to  the  bifurcating  limit 
cycles  as  the  Lewis  number  vanishes. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  shows  explicitly  the  effect  of  the  Lewis  number 
(a  measure  of  the  relative  importance  of  the  thermal  and  material 
time  scales)  on  the  dynamic  behaviour  of  a wide  variety  of  chemi- 
cally reacting  systems.  The  mathematical  structure,  which  is 
the  same  for  all  the  diverse  physical  systems  treated  (which 
include  stirred  tanks,  and  porous  and  nonporous  catalysts),  shows 
that  the  dynamic  state  tends  to  one  of  relaxation  oscillations  as 
the  Lewis  number  vanishes.  The  results  will  aid  the  understanding 
of  experimental  observations  in  many  chemical  systems. 


The  responsibility  for  the  wording  and  views  expressed  in  this 
descriptive  summary  lies  with  MRC,  and  not  with  the  authors  of 
this  report. 


THE  INFLUENCE  OF  THE  LEWIS  NUMBER  ON  THE 
DYNAMICS  OF  CHEMICALLY  REACTING  SYSTEMS 

4-  j. 

W.  Harmon  Ray'  and  Stuart  P.  Hastings* 

INTRODUCTION 

Chemically  reacting  systems  show  very  interesting  dynamics 
for  a wide  range  of  applications  [1-4].  The  influence  of  various 
parameters  on  this  behavior  has  been  studied  for  catalyst  parti- 
cles, continuous  stirred  tank  reactors  (CSTR's),  etc.  A particu- 
larly important  parameter  is  the  ratio  of  intrinsic  thermal  time 
constant  to  intrinsic  material  time  constant  (generically  referred 
to  as  the  "Lewis"  number)  for  these  reacting  systems.  It  is  well 
known  that  this  parameter  has  a very  striking  effect  on  the  system 
dynamics,  but  the  quantitative  effects  have  not  been  completely 
understood.  Previous  numerical  results  by  Hlavacek  et  al.  [5-7], 
Luss  and  Lee  [8,9],  and  others  [10-16],  showed  that  catalyst  par- 
ticles and  their  lumped  models  demonstrate  oscillations  for  suffi- 
ciently small  values  of  the  Lewis  number  in  some  cases.  Although 
these  results  have  been  particularly  well  treated  in  the  fine  sur- 
veys by  Aris  [15]  and  Sheintuch  and  Schmitz  [4],  the  form  of  the 
bifurcation  phenomena,  the  asymptotic  behavior  as  Le  0,  and 
the  quantitative  dependence  on  parameters  have  not  yet  been  dealt 
with . 

In  this  paper  we  endeavor  to  demonstrate  rather  explicitly 
the  influence  of  the  Lewis  number  for  several  major  classes  of 
reacting  systems  and  attempt  to  show  the  underlying  mathematical 
structure  which  transcends  a wide  range  of  specific  problems. 

First  we  develop  the  key  results  for  CSTR's  and  then  extend  these 
to  porous  and  nonporous  catalytic  surfaces.  One  particularly 
interesting  result  is  the  natural  appearance  of  relaxation  oscil- 
lations in  all  of  these  systems  as  the  Lewis  number  vanishes. 
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THE  CONTINUOUS  STIRRED  TANK  REACTOR 

The  dynamic  behavior  of  the  CSTR  has  been  analyzed  earlier 
by  bifurcation  theory  [17,18]  and  the  influence  of  reactor 
parameters  explicitly  determined  when  the  Lewis  number  was 
unity.  Here  we  wish  to  be  more  general  and  pose  the  CSTR 
modelling  equations  for  a first  order,  irreversible  exothermic 
reaction  as 


VM  — = F(CAf  " CA)  -Vkoexp{-E/RT}cA 
dt 


VT(PCp)T  = PCpF<Tf  " T)  + V ( - AH ) k Qexp { -E/RT } c^ 


-hA (T  - T ) (2) 

c 

The  quantity  VT^pCp^T  c°rresPonds  to  the  total  thermal  capacitance 

of  the  reactor  (perhaps  including  solid  packing,  reactor  wall, 

etc.)  while  VpC^  represents  the  thermal  capacity  of  the 

reacting  fluid  alone.  Similarly  V is  the  total  material 

M 

capacitance  of  the  reactor  while  V is  the  material  capacitance 
of  the  reactive  part  of  the  reactor.  As  an  example,  volume  V 
might  be  influenced  by  a catalyst  while  Vw  - V is  not.  A 
recent  experimental  study  [19]  illustrates  the  necessity  of 
this  more  general  formulation. 


By  defining  the  dimensionless  quantities 


- 3- 


one  obtains 


dx . 


— = -Xl  + Dad  - xx)  exp  u + , - 


= f.  ( x ,x„) 


(4) 


aX2  j X2 

Le  dt”  = “X2  + BDa(1  ~ V eXp  FTT , /Y| 


■B(X2  * X2c> 


(5) 


When  one  assumes  Le  = 1,  these  model  equations  reduce  to  those 
treated  in  [17,18]. 

To  begin  our  discussion  we  should  provide  a general 
definition  of  the  Lewis  number.  Earlier  definitions  (e.g.;  [20]) 

are  limited  in  scope  and  the  physical  interpretation  has  not 
always  been  made  clear.  Here  we  define  the  Lewis  number  as  the 
ratio  of  the  physical  transport  thermal  time  constant  to  the 
physical  transport  material  time  constant.  As  we  shall  show, 
this  allows  an  appropriate  Lewis  number  for  each  reacting 
system  with  all  playing  the  same  role  and  having  the  same 
physical  interpretation. 

For  the  CSTR,  the  Lewis  number  can  be  written 


V„  (pc  ) 

T K d T 
Le  = *■ 

VM  PCp 


v(pc  ) At 

T ^ p T 

F pc  At 
p 


Vc 


FAc 


(time  required  to  raise  the  tempera- 
I ture  of  the  entire  system  by  an 
(amount  At  through  inflow  and  out- 
jflow  of  fluid  at  flow  rate  F 

ftime  required  to  raise  the  concen-  ] 
Jtration  of  the  entire  material 
volume  an  amount  Ac  through  in- 
flow and  outflow  of  fluid  at  flow 
rate  F 

(6) 


I 


which  is,  of  course,  the  ratio  of  the  physical  transport  thermal 
and  material  time  constants. 

By  carrying  through  the  bifurcation  analysis  with  respect 
to  Le,  just  as  was  done  earlier  with  respect  to  Da  and  T 
[17,18]  and  augmenting  this  with  an  analysis  of  relaxation 
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phenomena,  one  can  see  the  entire  structure  of  the  dynamic 
behavior.  Obviously  steady  state  multiplicity  is  not  influenced 
by  Le,  so  that  we  shall  only  consider  stability  and  oscillatory 
behavior  here.  Let  us  divide  Eq.  (5)  through  by  Le  to  yield 


dx  (1  + B)  B [ X 0x 

= - x + — Da(l  - x ) exp  + = f (x  ,x  ) 

dt  Le  ^ Le  x 1 + x /y  Le 


This  puts  Eqs . (4,7)  in  the  same  form  as  the  equations  treated  in 

[17,18]  if  one  defines 


i + el  * b 

, B * , x , 


Thus  all  of  the  graphs,  calculations,  etc.  found  in  [17,18]  can 
be  used  directly. 

To  illustrate  these  results,  let  us  take  the  case  when 

x2c  = 0,  and  look  at  Figure  1 which  is  adapted  from  [17].  Some 

very  interesting  features  may  be  noted: 

(i)  Along  any  straight  line  through  the  origin,  the 

steady  state  behavior  x,  , x.  as  a function 

Is  2s 

of  Da  is  invariant.  For  straight  lines  below  the 
multiplicity  line  (M)  (curve  1) , the  steady  state 
behavior  is  in  the  form  of  Figure  2,  while  for 
straight  lines  above  line  M (curve  2),  the  steady 
state  behavior  takes  the  form  of  Figure  3. 

(ii)  The  range  of  dynamic  instability  behavior  of  the 

A 

steady  states  depends  on  B,  a as  shown  in  Figure  1. 
By  reducing  Le  while  keeping  all  other  parameters 

A 

fixed,  the  parameters  B,  a move  to  the  right  along 
these  straight  lines.  Curve  1 type  lines  would  begin 
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■oc  allow  limit  cycle  behavior  for  Le  values 
corresponding  to  point  A on  Figure  1,  where  the 
straight  line  enters  region  V.  Curve  2 type  lines 
would  begin  to  permit  limit  cycles  to  appear  at 
point  B where  the  straight  line  enters  region  III 
(or  VI  as  the  case  may  be) . 

(iii)  The  bifurcation  diagram  with  respect  to  Le  at 
constant  values  of  Da  should  have  one  of  the 
forms  shown  in  Figure  4.  Here  we  use  the  same 
notation  as  in  [17,181;  solid  lines  denote  stable 
steady  states  while  dashed  lines  are  unstable  steady 
states.  Solid  dots  represent  stable  limit  cycles 
while  open  dots  correspond  to  unstable  limit 
cycles.  The  particular  type  of  bifurcation  curve 
found  for  each  case  is  given  in  Table  1 and  depends 
strongly  on  the  sequence  of  regions  traversed  in 
Figure  1 as  well  as  the  value  of  Da  chosen.  Note 
that  stable  periodic  oscillations  seem  to  persist 
all  the  way  down  to  Le  = 0.  As  we  shall  show,  in 
some  cases  these  harmonic  oscillations  become 
asymptotic  to  relaxation  oscillations  as  Le  -*■  0. 
However,  no  relaxation  oscillations  are  possible 
for  the  case  when  y •*  00 . 

In  order  to  see  these  ideas  more  clearly,  let  us  linearize 
(4,7)  about  the  steady  state  to  obtain 

= A y 


Eqs  . 


dy 

dt 


(9) 


where 


V 


A 


3fl 

1 " 

! Hi 

3xl 

3X2 

3f2 

3f2 

3X1 

3X2 

x.  , x ^ 
Is  2s 


As  noted  in  [17,18],  the  bifurcation  points  to  periodic  solutions 
are  the  real  roots  of  tr  A = 0 which  also  satisfy  Det  A > 0 

r*S  'V-' 

and  we  can  examine  the  locus  of  these  points  as  Le  varies  by 
solving  the  cubic 


Here 


C 3 2 

tr  £ = 17  Cs  + a2s  + ais  + a0] 

R ^ 


= 0 


r = 


C = 


1 + 3 


r (1  + 3) 


Y2U-Xis)  (1  +x2s/Y)2 


Le 


a0  1 + 1 + 3J 


'4 

r 


r r r 


1 + 


Le 


1 + 3J 


a 2 = 


1_  + + i _ iX 

r 1 + e r 


J 


The  cubic  in  (10)  may  be  solved  numerically  to  determine  the 
bifurcation  points,  si,S2‘  The  crltlcal  value  of  Lewis  number, 
Le^,  above  which  there  are  no  oscillations,  is  the  solution  to 


(2a2  ~ 9a2ai  + 27a0)2 


[3a.  - (a  _ ) ] 


2,3 


0 


which  has  the  solution 


is  required  for  real  roots.  Thus  Le^  in  this  case  is  given  by 


Le 

c 


(1+6) 


In  general  for  values  of  Le  < Lec»  there  is  a range  of 
over  which  periodic  solutions  will  exist  provided  Det  A > 0.  Th 
actual  bifurcation  value  of  Le  depends  on  the  value  of  Da 
selected . 

The  locus  of  the  bifurcation  points  si'S2  are  expHcltly 

calculated  for  two  cases  in  Figures  2 and  3 and  these  serve  to 

illustrate  the  essential  points.  Note  that  for  parameters 

4 Y 

yielding  a unique  steady  state  always,  i.e.  r < ___  as  in 

y_4 

Figure  2 and  for  parameters  allowing  multiple  steady  states  (i.e 
r > as  in  Figure  3),  the  range  of  transient  instability 

s < x < s expands  as  Le  ■+  0.  This  may  be  seen  explicitly 
for  the  case  y ■+  where  s^  + 1/r  and  s 2 -*•  1 as  Le  + 0. 


Thus  for  a fixed  value  of  Da, 


limit  cycles  are  expected  to  persist  down  to  very  low  values  of 
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Le  as  shown  in  Figure  4.  Ir.  fact  for  the  case  of  finite  y, 
these  limit  cycles  become  asymptotic  to  relaxation  oscillations 
as  Le  vanishes.  This  may  happen  both  fo~  the  case  of  a unique 
unstable  steady  state  (Figure  5)  or  for  the  case  of  three 
unstable  steady  states  (Figure  6).  The  relaxation  trajectory 
begins  at  point  A (Figures  5,6),  follows  the  quasi-steady-state 
defined  by  = 0,  > 0 to  point  E where  it  must  jump  quickly 

to  C,  follows  x2  = 0,  < 0 again  to  point  D where  it  jumps 

quickly  to  A again.  Clearly  this  process  repeats,  giving  a 
periodic  solution  of  the  limiting  or  reduced  problem 


dt 

0 


flUl 


f2(Xl 


) 

) 


It  is  useful  to  discuss  the  properties  of  these  relaxation 
oscillations  (c'f . Stoker  [21]  for  a more  general  treatment). 

Firstly,  it  is  easy  to  show  that  relaxation  oscillations  do  not 
exist  for  the  case  y -*■  00  because  the  isocline  x^  = 0 does  not 
have  the  proper  form  (i.e.  three  branches  do  not  exist). 

Secondly,  for  relaxation  oscillations,  it  is 

necessary  that  the  steady  states  of  the  system  (corresponding  to 
intersections  of  x^  = 0,  &2  = 0)  fall  on  the  middle  branch  of 
x2  = 0.  It  is  easy  to  show  as  well  that  for  steady  states  on  either  the 
upper  or  lower  branch  of  x = 0,  tr  A < 0 always  and  no  bifurcating 
periodic  solutions  are  possible.  Thus  whenever  there  are  limit  cycle 
oscillations  around  a unique  steady  state  or  around  three  unstable 
steady  states,  these  may  become  asymptotic  to  relaxation  oscilla- 
tions of  the  same  type. 


i 


Figure  2 Steady  state  behavior  and  locus  of  bifurcation 
points  as  Le  changes;  y = 20,  r = 4.05,  B = 


3- 


pes  of  Bifurcation  Behavior  for  Le 


Types  of  Bifurcation 
(Figure  4) 


(i)  Type  (b)  for  left  hand  bifurca- 
tion point  in  region  Va 


(ii)  Type  (a)  otherwise 


2.  I -*•  VI  -*■  Illb  -*•  IV 

(i) 

Type  (d)  for  bifurcation  in  the 
region  of  multiplicity 

(ii) 

Type  (a)  for  the  upper  branch 
in  region  of  uniqueness 

(iii) 

Type  (b)  for  the  lower  branch 
in  region  of  uniqueness 

(i) 

Type  (c)  or  type  (d)  for 

bifurcation  in  the  region  of 
multiplicity  (depending  on 
the  particular  value  of  Da) 

(ii) 

Type  (a)  for  the  upper  branch 
in  the  region  of  uniqueness 

(iii)  Type  (b)  for  the  lower  branch 
in  the  region  of  uniqueness. 
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As  an  illustration  of  these  points,  consider  Figures  5 
and  6 where  computed  stable  limit  cycles  are  shown  as  a 
function  of  Lewis  number.  The  case  of  a unique  steady  state 
corresponding  to  curve  1 in  Figure  1 and  Da  = 0.15,  is 
shown  in  Figure  5.  For  Le  = 1.0,  the  steady  state  is 
stable  and  surrounded  by  both  stable  and  unstable  limit 
cycles  (to  avoid  cluttering  the  drawing,  only  the  stable  one 
is  shown).  As  Le  decreases  to  0.5,  the  steady  state  becomes 
unstable  and  the  limit  cycle  becomes  very  large.  Further 
reduction  in  Le  causes  the  limit  cycle  to  approach  the 
relaxation  oscillation.  Clearly  this  example  is  a case 
of  type  (b)  bifurcation  in  Figure  4. 

For  the  case  of  three  steady  states,  the  dynamic  behavior 
as  Le  varies  along  curve  2 in  Figure  1 was  computed  for 
Da  = 0.090.  The  bifurcation  behavior  for  this  case  is  of 
type  (c)  in  Figure  4.  From  about  Le  = 2.25  down  to 
Le  = 2.16,  an  unstable  limit  cycle  was  found  surrounding  a 
stable  upper  steady  state.  From  Le  = 2.16  to  about 
Le  = 0.33,  only  the  lower  steady  state  is  stable  with  no 
limit  cycles  present.  At  about  Le  = 0.33  an  unstable 
limit  cycle  appears  around  the  lower  steady  state  and  a 
stable  limit  cycle  appears,  encircling  all  three  steady 
states.  For  Le  < 0.21,  only  the  large  stable  limit  cycle 
remains  and  it  grows  asymptotic  to  the  relaxation  oscillation 
as  Le  -*  0 . 
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NOMPOROUS  CATALYTIC  SURFACES 

For  nonporous  catalytic  material s such  as  catalytic  wires 
or  gauzes,  the  dynamic  equations  can  take  the  general  form  of 
Eqs.  (4,7)  where  depending  on  the  specific  model  a s sumptions , 
the  parameters  Da,  B,  y,  6,  Le  can  have  several  meanings  (22). 
To  illustrate  a specific  case,  consider  the  model  of  Ervin  and 
Luss  [23]  for  an  infinite  length, cylindrical , catalytic  wire  upon 
which  oxidation  reactions  are  carried  out.  This  model  assumes 
that  mass  transport  to  the  wire  surface  is  controlled  by 
irreversible  adsorption  of  reactant,  and  that  the  reaction 
products  have  negligible  surface  coverage.  In  primitive  form, 
the  model  equations  are 


T - Vg^T  " V "k0  6 


-E/RT 


ir(d  ) dT 

p C — - * TTd  h (T  - T ) + ( -Ah ) 7r d k . e 

. wg  s wO 

4 d t 


-E/RT 


where  9^  and  0T  represent  the  surface  coverages  of  reactant 

and  total  surface  respectively,  T , c are  the  gas  phase 

g 9 

temperature  and  reactant  concentration,  and  k^  represents  a 
temperature  independent  adsorption  rate  constant.  Further,  h 
is  a heat  transfer  coefficient,  Tg  the  wire  temperature,  and 
d^  the  wire  diameter.  By  defining  the  dimensionless  quantities: 


8_  - 6. 

T A 


T - T k e 

— 2-  (E/RT  ) , Da  « — 

T 9 k,c 

g i g 


t = t k c , y = E/RT  , B 
i.  9 9 


(-AH)klcq6TY 


, Le 


do  c k c 
w^  s ps  1 c 


e = o 
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then  tqs . (17,18)  take  the  form  of  Eqs  . (4,7).  Here  our 

definition  or  the  Lewis  number  for  the  catalytic  wire  is  totally 


consistent  with  our  general  definition  presented  above,  for 

( IT  / J \ 2 1 


d p C k c 

w s ps  1 < 


"(dw> 


p c At/it d hAT 

s ps  w j 


(0_,/k  c 8„,) 
T 1 g T 


time  required  to  raise  the  wire  temperature] 
an  amount  AT  through  heat  s 
.transfer  from  the  gas J 

time  required  to  increase  the  surface  concentra 


by  an  amount 


due  to  adsorption  on  the  surfa 


tionl 
ace  / 


which  is,  of  course,  the  proper  ratio  of  physical  transport  time 
constants . 

Having  put  (17,18)  into  the  equivalent  CSTR  form  (4,7), 
all  of  the  analysis  of  the  previous  section  immediately  applies. 
Thus  one  may  readily  determine  the  influence  of  the  Lewis  number 
(as  well  as  other  parameters)  on  the  dynamic  behavior.  From 
a practical  point  of  view  one  may  calculate  the  critical 
Lewis  number,  Lec  and  compare  with  the  actual  value  to 
determine  whether  or  not  periodic  solutions  or  relaxation 
type  oscillations  are  to  be  expected.  As  noted  in  [22],  for 
the  most  common  models  of  catalytic  wires,  such  oscillations 
are  not  predicted  by  the  equations.  However,  there  are  models 
for  which  Le  is  in  the  proper  range  to  allow  oscillations 
for  nonporous  catalytic  surfaces  [24]. 


POROUS  CATALYTIC  MATERIALS 

Porous  catalytic  materials  such  as  catalyst  particles, 
catalyst  impregnated  gels,  etc.  are  known  to  have  exotic 
dynamics  which  are  influenced  by  the  Lewis  number.  To 


4 
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illustrate,  let  us  consider  the  case  of  a single  irreversible 
first  order  reaction,  for  which  the  basic  equations  take  the 
form  [15] 


, -E/RT 

_ko  e c 


PC. 


£ 


and 


D c 

& 3 1 

a 

z 

3 c 

3 1 ' 

za  3 z i 

3 z 

3t  k 3 

* 

a 

2 

3 T 

at' 

a ~ 

Z d Z 1 

i 

3 z 

9 c 

d z 

= k (c 
m g 

c ) 
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where  c is  the  reactant  concentration,  T the  temperature, 

the  mass  diffusivity,  k the  thermal  conductivity,  k and  h 

m 

external  heat  and  mass  transfer  coefficients,  and  a is 
geometric  parameter  having  the  following  values:  a = 0 (slab), 
a = 1 (cylinder),  a = 2 (sphere).  By  defining  the  dimensionless 
quantities 


c - c 
x.  = -2 


T - T 
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RT 
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, r = z/i 
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, Le  = 


E/RT 


(25) 


kT 


hX. 


Nu  = 


& k 

one  obtains  the  dimensionless  equations 
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dx 

3x  1 3 r*  — * — _ [ x 

— - t “ — - + ♦ (i  - v exp  — — rr  gi(xi'x2)  (26) 

3i  r 3r  | 1 + x2/y  j 

3 x ] 

3x2  1 nr#3ri  2 x2 

Le  = — + B<f>  (1  - x ) exp  = g (x  , x ) (27) 

3t  ra  3r  1 + x2/\ 


+ Shx, 


+ Nux,  = 


r = 1 


3xx  3x2 
3 r 3 r 


r = 0 


Notice  that  (28,29)  are  Robin  type  boundary  conditions;  however 
as  Hu,  Sh  they  reduce  to  Dirichlet  conditions  of  the  form 


*2  = 0 ; 


Consistent  with  the  earlier  definition,  the  Lewis  number  for 
this  catalyst  particle  is  the  ratio  of  appropriate  transport  time 
constants ; i . e . , 


*C.E  . 


v pc  At 

p p 

A kAT /l 


V Ac 
P 

A & Ac/Z . 
P ' 


time  required  to  raise  the  catalyst 
particle  temperature  an  amount  At  due  to 
heat  conduction  across  the  external  surfacej 

time  required  to  raise  the  catalyst  parti- 
cle concentration  an  amount  Ac  due  to  L 

diffusive  mass  transport  across  the 
external  surface 


By  the  use  of  any  number  of  lumping  procedures  [15],  Eqs . 
(26-30)  may  be  reduced  to  an  equivalent  set  of  CSTR  equations  of 
the  form  of  Eqs.  (4,7).  Thus  the  effect  of  Lewis  number  may 
be  seen  by  using  the  CSTR  results  developed  above  together  with 


one  of  these  lumping  procedures. 


To  provide  an  example,  consider  the  lumping  procedure  of 


Llavacek  and  coworkers  [5-7]  for  the  case  of  a slab  geometry 
(a  = 0 in  (26,27)).  Their  equivalent  stirred  tank  equations 
take  the  forn 


dx 


dx 


— = -(X^2^  + ~ ) exp| 


i1  + x2/T. 


^X2  2 2 
Le  = -(U,)  x + B4>  (1  - x ) exp < 

dx 


1 + x2/y 


(33) 


(34) 


where  X^,  y are  first  positive  roots  of  the  transcendental 
equations 


tan 

X = 

- X/Nu 

(35) 

tan 

V = 

-y/sh 

(36) 

Note  that  as 

Nu, 

Sh  -»•  », 

and  the  Dirichlet  boundary 

conditions 

apply,  X2  = 

U1  " 

IT.  By 

defining 

t = (X^  x , Da  = 


(XL) 


T ' 1 + 6 ■ 'V2/(V2  - x2c  ■ 0 


Equations  (33,34)  are  exactly  in  the  form  of  the  CSTR  equation 
(4,7).  Thus  all  of  the  earlier  CSTR  analysis  and  resulting 
structure  also  applies  to  the  catalyst  particle  problem. 

In  spite  of  numerical  calculations  showing  reasonable 
agreement  between  the  behavior  of  the  distributed  and  approximate 
lumped  systems,  the  formal,  nonrigorous  nature  of  these  lumping 
approximations  to  distributed  reacting  systems  may  cause  concern. 
Thus  it  is  useful  to  consider  more  direct  evidence  that  suggests 
we  have  the  correct  structure.  Firstly,  a number  of  detailed 
numerical  case  studies  of  the  distributed  catalyst  particle 
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[6,8-14,16]  with  varying  Lewis  number  show  the  type  of  behavior 
predicted  from  the  lumping  analysis.  In  particular,  the  oscilla- 
tions observed  become  stiffer  with  decreasing  Le  so  that 
distributed  relaxation  oscillations  seem  to  be  the  plausible 
asymptote.  The  general  nature  of  these  distributed  relaxation 
oscillations  can  be  seen  with  the  aid  of  Figure  5 and  Eqs.  (26-30)  . 
When  Le  -►  0,  then  92^Xl'X2^  = ° and  together  with 

appropriate  boundary  conditions  define  the  behavior  of  the 
system.  For  a given  composition  profile,  x^(r),  the  equation 
g^(x^,x2)  = 0 and  associated  boundary  conditions  could  have 
three  quasi-steady-state  solutions  x^(r),  x”(r) , x^(r)  as 
shown  in  Figure  7.  This  is  similar  to  the  situation  for  the 
lumped  system  shown  in  Figures  5 or  6 where  for  any  value  of 


x^  between  points  A and  B,  there  are  three  quasi-steady-state 

;lutions  x2>  It  is  conjectured  that  by  analogy  with  the  lumped 
3x1 (r , t)  3x  (r , t) 

case,  <0  on  the  lower  branch  X2£^r^'  and  9t > 0 

u 

on  the  upper  branch  x (r),  so  that  a distributed  relaxation 

£ 

oscillation  jumping  between  a sequence  of  lower  solutions  x2(r) 
and  a sequence  of  upper  solutions  x2(r)  should  occur. 

Further  evidence  of  this  behavior  may  be  seen  if  we  let 


the  composition  profile  x^(r)  be  the  steady  state  solution 

corresponding  to  a steady  state  value  of  x2ss  (r)  on  the  middle 

$L  m 

branch  such  as  in  Figures  5 and  6.  Recalling  that  x (r)  < x_,(r) 
u 3*i<r,t) 

= x2ss  (r)  < x2(r),  we  can  see  that  ^ in  Eqn.  (26)  will 

£ 

have  the  signs  noted  above  if  x (r)  is  replaced  by  x (r) 

4 S S 2. 

(lower  branch)  or  x2(r)  (upper  branch).  This  is  due  to  the 
monotone  increasing  nature  of  the  term  exp  ( ^ ^ ^ ^'/y  } Ecjn* 


" 
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One  should  note  that  the  observed  oscillations  of  Lee  and 
Luss  [8]  for  Le  = 0.1  are  rather  stiff;  however,  their  spatial 
character,  showing  a transient  hot  spot  close  to  the  external 
surface,  cannot  be  due  to  purely  relaxation  phenomena.  This 
illustrates  the  point  that  relaxation  oscillations  must  be 
considered  no  more  than  an  asymptote  to  the  bifurcating  limit 
cycle  oscillations  found  at  finite  values  of  Lewis  number. 

CONCLUDING  REMARKS 

It  has  been  suggested  that  the  Lewis  numbers  necessary 
to  allow  oscillations  in  catalytic  systems  are  much  too  small 
to  be  realistic.  However,  as  can  be  seen  from  recent  work 
(e.g.,  [4,24]),  more  detailed  models  indicate  that  such  small 

Lewis  numbers  may  be  realized  in  practice  and  consideration  of 
Lewis  numbers  close  to  zero  is  a fruitful  exercise. 

A conjecture  put  forth  by  a number  of  workers  is  that 
for  Lewis  numbers  smaller  than  some  bifurcation  value,  the 
steady  state  in  question  becomes  unstable  and  oscillations 
arise,  while  for  values  greater  than  this  bifurcation  value  the 
steady  state  is  stable  and  no  oscillations  arise.  The  present 
analysis  shows  that  this  conjecture  is,  in  principle,  correct; 
however  as  seen  in  Figure  4,  things  are  a little  more 
complicated.  For  a unique  steady  state,  bifurcation  can  be 
to  the  right  as  in  Figure  4(b)  and  stable  oscillations  observed 
even  when  the  steady  state  is  stable.  Furthermore,  for  the  case 
of  multiple  steady  states  as  in  Figure  4 (c) , (d) , there  are 
different  bifurcation  values  of  the  Lewis  number  for  the  upper 
ana  the  lower  branches. 

Although  a rather  wide  range  of  reacting  systems  have  been 
considered  here,  it  has  been  possible  to  put  them  all  within 
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unified  mathematical  structure.  This  has  the  advantage  that  in 
analyzing  new  problems  of  stability  and  oscillations,  these 
results  may  serve  as  an  aid  to  physical  intuition  in  order  to 
evaluate  the  dynamic  phenomena  to  be  expected. 
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